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The structural and electronic properties of a single pentacene molecule and a pentacene molecular 
crystal, an organic semiconductor, are examined by a first-principles method based on the generalized 
gradient approximation of density functional theory. Calculations were carried out for a triclinic unit 
cell containing two pentacene molecules. The bandwidths of the valence and conduction bands which 
determine the charge migration mechanism are found to depend strongly on the crystallographic 
direction. Along the triclinic reciprocal lattice vectors A and B which are orientated approximately 
perpendicular to the molecular axes the maximal valence (conduction) band width amounts to only 
75 (59) meV, even smaller values are obtained for the C direction parallel to molecular axes even 
less. Along the stacking directions A+B and A-B, however, the maximal valence (conduction) 
band width is found to reach 145 (260) meV. The value for the conduction band width is larger 
than estimates for the polaron binding energy but significantly smaller than recent results obtained 
by semiempirical methods. The single molecule has a HOMO-LUMO gap of about 1.1 eV as deduced 
from the Kohn-Sham eigenvalue differences. When using the self-consistent field method, which is 
expected to yield more reliable results, a value of 1.64 eV is obtained. The theoretical value for the 
band gap in the molecular solid amounts to 1.0 eV at the F-point. 



I. INTRODUCTION 

Organic semiconductors based on pentacene {C22H14) 
or other aromatic hydrocarbon molecules have recently 
attracted an enormous interest regarding their use for 
molecular electronics Q]. A pentacene field effect tran- 
sistor showing a temperature-independent high-mobility 
(/i ^ l.icni^ /V s) has been fabricated[Q. It has to be 
noted, however, that in the latter case polycrystalline 
material was used; considerable higher mobilities can be 
expected for single crystalline material. Although mo- 
bilities for single crystalline pentacene have not yet been 
measured, recent work has demonstrated that epitaxial 
growth of pentacene on solid substrates may be possible 
thus putting the realization of devices with an ac- 
tive layer of single-crystalline pentacene into sight. De- 
spite this pronounced interest, however, there are still 
many open questions even on rather basis properties 
on molecular semiconductors like pentacene, in partic- 
ular as far as the electronic structure is concerned. The 
conventional view is that in molecular crystals of aro- 
matic molecules like pentacene hole transport is domi- 
nant and that the valence band width is of order 100 
meV - much smaller than the experimental 500 meV 
estimate for the conducting charge transfer salt TTF- 
TCNQ (tetrathiafulvalene-tetracyanoquino dimethane). 
The bandwidth is further important for understanding 
the charge transport mechanism. According to Holstein's 
polaron modelQ, the small polaron limit (size of lat- 
tice constant) is reached when the electronic bandwidth 
is small compared to the polaron binding energy. The 



binding energy is experimentally estimated to be of or- 
der ~ 200 meV[^. On the theoretical side, however, 
there have only been very few reports. Using semiem- 
pirical methods, results ranging between ~ QQQmeV ob- 
tained using the cluster approach with a semiempirical 
Hartree-Fock INDO (intermediate neglect of differential 
overlap) and ~ 12QmeV obtained from an extended- 
Hiickel-type (EHT)-calculation[|| were reported. 

In this paper, we address the above mentioned issues 
by using a first-principles algorithm In particular, we 
illustrate how the bandwidths of the VB and CB depend 
on the orientation of the herringbone structured crystal. 
We find that although there is only weak dispersion of 
the bands along the reciprocal lattice vectors A, B, and 
especially C (along the molecules), along the stacking 
directions A+B and A-B the dispersion reaches signifi- 
cantly larger values. For these directions, the bandwidth 
of the CB and the next higher band CB+1 amounts to 
149 and 260 meV. For the A-B direction the VB and the 
next lower band VB-1 don't become quite as large, the 
corresponding values amount to 145 and 131 nieV. 

Hence, the band widths in the stacking directions can 
become larger than estimates for the polaron binding en- 
ergy. This should make band-like transport at low tem- 
peratures possible. To the best of our knowledge this 
is the first fully ah initio bandstructure calculation and 
theoretical band gap estimate in the literature. 

In section ^ we discuss briefiy the method of calcu- 
lation. It will be followed by a section describing the 
models in section pl[ Results of the single molecule and 
the molecular solid and discussions will be presented in 
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FIG. 1: Atomic bond-lengths a) of a fully relaxed single pen- 
tacene molecule (in xy-plane) and b) of the constrained re- 
laxed pentacene crystal. One can clearly see the alternating 
double (shorter) and single (longer) bonds. The experimen- 
tal values are shown in brackets (...). The missing bond- 
lengths follow from symmetry. 



section [V. Finally, in section 0, a summary will be 
given. 

II. METHODS OF CALCULATION 

We used one of the popular density functional the- 
ory (DFT) algorithms with localized orbitals as the ba- 
sis functions, SIESTAQ. It uses TrouUier-Martins 
norm-conserved pseudopotentials Jill in the Kleinman- 
Bylander separable form[^. The basis set is made of 
pseudo atomic orbitals (PAO) of multiple-zeta form in- 
cluding optionally polarization orbitals. The first-zeta 
orbitals are obtained by the method of Sankey and 
Niklewski[p^, while the second-zeta orbitals are con- 
structed in the split-valency philosophy well known from 
quantum chemistry With this basis set, SIESTA 
calculates the self-consistent potential on a grid in real 
space. The fineness of this grid is determined in terms of 
an energy cutoff Ec in analogy to the energy cutoff when 
the basis set involves plane waves. In the present calcu- 
lations, we used Ec to be 80 Ry and a double-zeta plus 
polarization orbitals (DZP) basis set. For the exchange- 
correlation energy functional, the generalized gradient 
approximation (GGA) in the version of Perdew-Burke- 



III. MODELS 

A. Single molecule 

For the case of simulating a single molecule, we used 
a large cubic supercell {50A x 50A x 50A) and placed 
the molecule at the center. These dimensions are much 
larger then the pentacene dimensions (14.21A x b.QAA 
in x-y plane) and results in no interactions between the 
unit cells. The atoms in the molecule were relaxed until 
the magnitude of the force on each atom is less than the 
tolerance OMeV/A. 



B. Molecular solid 



For the molecular solid, we used a triclinic unit cell 
containing two non-equivalent molecules. The unit cell 
data with the Bravais lattice vectors a, b (perpendicular 
molecular axes) , and c (along molecular axes) and atomic 
positions were taken from a previous x-ray measurement 
on a trichlorobenzene solution-gown single crystals by R. 
Campbell and co-workers (1961) The lattice vec- 
tors are a = 7.93, b = 6.14, c = 16.03A, a = 101.9, 
f3 = 112.6, 7 = 5.8°. There are other structure determi- 
nations, e.g. by D. Holmes et al. (1999)0. Although 
the crystallization solvent was the same in both cases, the 
unit cell parameters disagree slightly. Their structure is 
similar to crystals obtained from the vapor-phase depo- 
sition technique which results in slightly denser molecu- 
lar packing (unit cell volume reduced by 3%)|l^. Poly- 
morphism is common in molecular crystals and depends 
on the applied preparation techniques and conditions Jl8[. 
The main structural features however are the same and 
we restrict ourselves to the analysis of pentacene crystal 
described by Campbell's data. 

Since the GGA functional does not include van der 
Waals attraction! 19 , 20) , we relaxed the atoms except 3 in 
each molecule in order to fix the planes of the molecules. 
The force tolerance is again O.OieV/ A. The k-point sam- 
pling is done just using the F-point. This is expected to 
be sufficient due to the large unit cell. 



In k-space, we denote the reciprocal lattice vectors A, 
B, and C. We note however that the unit cell is triclinic 
and not tetragonal, although the unit cell angles are close 
to 90°. Hence by convention, e.g. A is parallel to b x c 
but not to a. 
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FIG. 2: Single particle wavefunctions of isolated pentacene 
molecule. The dark color represents a positive sign, the light 
color a negative sign, (a) and (b) show the HOMO, (c) and 
(d) the LUMO. (a) isosurface at ±0.05/1^/^; (b) LCAO co- 
efficients: small circles O.Of - 0.09, large circles 0.12 - 0.3, the 
short lines show regions of high gradients; (c) isosurface at 
±0.03/1^/^, (d) similar to (b). Wavefunctions plotted with 
the program gOpenMol|24]. 



IV. RESULTS AND DISCUSSIONS 

A. Single molecule 

In figure ^ a) we show the relaxed single molecule. 
We label the atoms and give the bond- lengths. The z- 
coordinates for the relaxed atoms are less than 0.67mA 
about z = showing that the molecule is essentially pla- 
nar, in agreement with the experimental findings [ITj. 
The molecule is symmetric with respect to the C atoms 
6 and 17 and the midpoints of the bonds 1-22 and 11-12. 



The bond- length in the y-direction is typically 1.46A, 
except at the ends where it is 1.43 A. We will compare 
the C-C and C-H bond-len gths to experimental crys- 
tallographic data in section 



IVB 



In general they are 
similar to those in graphite (1.42A). Our results show 
also smaller bond-lengths for the C atoms near the ends 
(bonds 1-2, 10-11, 12-13, and 21-22) to their respective 
neighbors along the x-direction. These bond-lengths are 
1.381. All the C-H bonds are l.lOl long. 

The energy gap A separating the HOMO and the 
LUMO states is 1.1 eV as deduced from Kohn-Sham 
eigenvalue differences. This value is about 40% smaller 
than the experimental value, 1.82 eV, reported from el- 
lipsometric spectra of thin pentacene films . This un- 
derestimation of gaps by this magnitude is a common 
deficiency of the GGA functional |^^. Since DFT is 
a groundstate theory, better results can usually be ob- 
tained from total energy differences (self-consistent field 
method) H 



A^Eo{N)-Eo{N + l), 



(1) 



where Eo{N) is the groundstate of the neutral N (even) 
electron system and Eo{N + 1) is the groundstate with 
one extra electron added in the same geometry as the 
neutral system. This gives a quantitatively better gap 
estimate of 1.64 eV. 

The corresponding wave functions for these two states 
are shown in figure ^ Surprisingly, they are different 
from a previous calculation by Strohmaier et al. ||2^ us- 
ing semiempirical MNDO (modified neglect of diatomic 
overlap). Their HOMO resembles our LUMO, while their 
LUMO has the same symmetries as our HOMO although 
contributions from carbon atoms near the end of the pen- 
tacene molecule seem different. We do not know the rea- 
son for the discrepancy. However, we checked our result 
further by comparing it to an elementary Hiickel-type 
calculation only including the 22 orbitals each con- 
tributing one electron. The signs of the wavefunction are 
the same as the ones obtained from the more elaborate 
SIESTA calculation. 

Fig. I (a) shows the HOMO state. The main contri- 
bution to this wave function comes from the Pz orbital 
of the C atoms. In Fig. || (b), we show the coefficients 
of the wave function expanded onto the first-zeta basis 
functions which resemble most closely the atomic orbitals 
with pz symmetry. The positive values are indicated by 
the filled circles and the negative values are shown by the 
open circles. The large circles correspond to a magnitude 
between 0.12 and 0.3, while the small circles represent a 
value between 0.01 and 0.09. Both figures clearly demon- 
strate that the C atoms contributing to the HOMO state 
are in a second neighbor (two large circles with a small 
circle in between) configuration near the center but in 
the nearest neighbor configuration at the two ends of the 
molecule. It is important to note that the relative phase 
of this molecular wave function varies by 180° in the near- 
est neighbor configurations at the ends of the molecule. 
In Fig. H (c), we plot the wave function of the LUMO 
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FIG. 3: Herringbone packing of two inequivalent pentacene 
molecules in the unit cell of Campbell's model. The molecules 
are relaxed although three atoms in each molecule are kept 
fixed. The Bravais lattice vectors a, b, and c are assigned, as 
well as the vectors a+b and a-b. Figure prepared with the 
program MOLMOL|||]. 

state. Similar to the HOMO state, the dominant eon- 
tributing orbitals are also the Pz orbitals. As shown in 
Fig. ^ (d), similar second neighbor configuration near 
the center as in HOMO state is observed but there are 
more nearest neighbor configurations at the ends of the 
molecule. This difference at the ends of the molecule ex- 
plains why the HOMO state has lower energy than the 
LUMO state because the latter state has a larger kinetic 
energy due to rapid (180°) phase changes. 

B. Molecular solid 

In figure |^ we show the structure for the molecular pen- 
tacene solid according to the results of Campbell. The 
triclinic unit cell for each case is depicted by thin solid 
lines. There are two molecules per unit cell labeled I and 
II following Campbell's notation which are not parallel 
to each other, but are packed in a herringbone fashion. 
The comparison of our atomic bond- lengths from the con- 
strained relaxation with the experimental ones (in brack- 
ets) are shown in figure ^ b) . The centers of symmetry 
of molecule I (0,0,0) and molecule II (1/2 ,1/2 ,0) in the 
space group PI are marked by arrows. The two redun- 
dant bond-lengths at the right end of molecule I illus- 
trate this symmetry. The bond-lengths are qualitatively 
correct and show the proper trends (alternating single- 
double bonds), but are in general slightly larger than the 
experimental ones as typical for the DFT method. Addi- 
tionally, the experimental bond-lengths are not symmet- 
ric with respect to the long molecular axes as opposed 



to the DFT result where the bonds are more symmetric 
similar to the single molecule (at least to this precision). 
This reflects that noncovalent bonding by nonlocal inter- 
actions is not treated properly. This deficiency however 
concerns mostly the total energy (and structure) , less the 
single particle wavefunctions (and bandstructure) . 

The band structure along the high symmetry axes for 
energies ranging from -8 to -3.8 eV are depicted in figures 
^ and H The bands always come in sets of two, due 
to the near degeneracy of the states of two pentacene 
molecules per unit cell. There are two important pieces 
of information one can read off the band structure plots. 

First, we want to estimate the transfer integrals t be- 
tween the two molecules in a unit cell and restrict our- 
selves to the couplings between the two HOMOs and be- 
tween the two LUMOs. These are most important for the 
hole and electron transport, respectively. Since there is 
no momentum transfer involved for the intracell coupling, 
we can restrict ourselves to the F-point. At this point 
the bonding/antibonding splitting between two (identi- 
cal) HOMO (LUMO) states is 2tH(L)- This gives tn = 85 
meV for the HOMO coupling and fz, = 15 meV for the 
LUMO coupling. Since the two molecules are in princi- 
ple inequivalent where the level energy difference between 
the two HOMO (LUMO) states is AEh{l), the true cou- 
pling i'ff(i) follows from 2t'^^^^ = ^{2tH(L)? - 

This equation of a coupled inhomogeneous two level sys- 
tem was stated by Cornil et (their Eq. (1)). In 
their work an offset energy AEu(^]^-^ as large as 61 meV 
(70 meV) was quoted. While Is.Eh = 61 meV leads to 
t'jj = 80 meV, Ai?L — 70 meV gives a negative discrim- 
inant, i.e. our and Cornil's Ai?i are not consistent. 
This at least illustrates that there cannot be a large offset 
in our case. 

Second, the dispersion and bandwidth reflect the cou- 
pling between the unit cells and determines the crystal 
properties. The larger the bandwidth the more delocal- 
ized the electronic states are at finite temperature and 
the more one has transport by a band mechanism. If 
the bandwidth becomes small and comparable to the po- 
laron binding energy, then excess charges are being self- 
trapped and need thermal activated to migrate by a hop- 
ping mechanism. 

The bandwidths are shown in table |. At least along 
the triclinic reciprocal lattice vectors A, B, and C they 
are all smaller than 73 meV and hence clearly smaller 
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TABLE I: Bandwidths of valence band (VB), the next lower 
band (VB-1), the conduction band (CB), and the next higher 
band (CB-H) in meV. 
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FIG. 4: Bandstructure along reciprocal lattice vectors A, B, 
and C. The LUMOs of the two molecules per unit cell couple 
less than the HOMOs, also the conduction band shows hardly 
any dispersion. 
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FIG. 5: Bandstructure along the stacking directions A+B 
and A-B. In these direction the bandwidths exceeds the po- 
laron binding energy and band-like transport should be pos- 
sible. 



than estimates of polaron binding energies (~ 200meV 
1^). This is likely to lead to a charge migration mecha- 
nism of some sort of hopping in these directions. Either 
the small polaron or the multiple trapping and release 
model Q have been suggested. Our values for the band 
widths along A, B, and C are in semi-quantitative agree- 
ment with previous semiempirical ETH calculations by 
Haddon and co-workers [|| . 

It was suggested by Cornil and co-workers using ex- 
trapolation of results from finite clusters that the band- 
width is much larger 600 meV) in the stacking di- 
rections connecting the two inequivalent molecules. In 



their Fig. 1, di connects the two molecules within same 
unit cell and d2 between neighboring unit cells. Here we 
calculated the bandstructure from first-principles along 
the high symmetry directions A-B and A+B resembling 
Cornil's di and d2 directions. We obtain that the width 
of the CB and in particular the next higher band CB+1 
increases drastically to 149 and 26QmeV, respectively. 
On the other hand, we only observe an increase of the 
width of the VB and next lower band VB-1 in the A- 
B direction. At a first glance, they seem to be 211 and 
247 meV. However, there is a band crossing for the two 
highest occupied bands. Since the triclinic unit cell has 
no symmetries (except the identity), the crossing bands 
mix and open a 10 meV gap (not resolved in plot). Due 
to the avoided crossing the widths of VB-1 and VB are 
only 131 and 145 meV, respectively. 

If one assumes a cross-over between polaron-like and 
band-like transport at a bandwidth of 200 meV, band- 
like electron transport should be possible in the stack- 
ing directions. However, we do not observe bandwidths 
as large as those reported in Cornil's work Q. When 
measuring the conductivity of a pentacene crystal or ex- 
tracting the mobility from transistor characteristics, one 
should obtain thermally activated charge hopping as well 
as temperature-independent band transport behavior de- 
pending on the crystal orientation. Exactly this was ob- 
served by S. F. Nelson and co-workers (see their Fig. 3) 

i- . 

The question whether the width of the CB is smaller or 
equal to the VB cannot be answered completely. Let us 
look at figure 1^ first. The LUMO splitting and dispersion 
of the CB is clearly smaller than the HOMO splitting 
and dispersion of the VB. This in line with conventional 
expectations. On the other hand, figure || reveals that the 
bandwidths of the CB along A+B can indeed become 
larger than the VB bandwidth. In particular, the widths 
of CB-f 1 is drastically enlarged for A+B and A-B. 

We further confirm by our method that there is hardly 
any dispersion in the C direction approximately along the 
pentacene molecules. As already pointed out in Ref. 0| 
this leads to a quasi 2-dimensional character for charge 
transport. 

Finally, the fundamental band gap is 0.97 eV mea- 
sured at the F-point. This is slightly smaller than the 
Kohn-Sham HOMO-LUMO gap for the single pentacene 
molecule, 1.10 eV, due to band offsets from Brillouin zone 
folding and bonding/antibonding splitting. As demon- 
strated in the case of a single pentacene molecule the 
gap is drastically underestimated. 



V. SUMMARY 

In conclusion, we used an ab-initio approach to de- 
termine the electronic properties and in particular the 
band-structure of a molecular solid of pentacene. A ge- 
ometry optimization of the single molecule yields an es- 
sentially planar molecule in agreement with experiment. 
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The bond-lengths along the x-direction (the direction 
of the molecule) do show 0.06A difference between the 
bonds at the ends and the middle of the molecule. Both 
the HOMO and LUMO states originate from the pz or- 
bitals on second neighbors C atoms in the center and 
on nearest neighbors C atoms at the two ends of the 
molecule. The relative phases and the extent of overlap- 
ping between the neighboring pz orbitals differ for the 
two states. 

For the solid pentacene, calculations we carried out 
using the first-principles tight-binding code SIESTA 
for a pentacene molecular crystal with the experimen- 
tally determined herringbone structure together with in- 
tramolecular distances resulting from a constrained ge- 
omtry optimization. The resulting bond-lengths are in 
good agreement with experiment. The solid is predicted 
to be a large band-gap (> l.QeV) semiconductor with 
an maximal bandwidth for electron transport of about 
260 meV and a maximal bandwidth for hole transport 
of only 145 meV. It is found that the widths of the elec- 
tronic bands depends strongly on the crystallographic di- 
rection. Along the triclinic reciprocal lattice vectors the 



bandwidths are generally smaller than estimates based 
on the small polaron binding energy whereas. Along 
the stacking directions a significantly larger width is ob- 
served. The present maximum bandwidths are much 
smaller than previous results obtained using the semiem- 
pirical INDO approach. On the basis of the present the- 
oretical results a band-like transport of charges in high- 
quality pentacene single crystal with very high mobilities 
should be possible at low temperatures, as previously ob- 
served experimentally for napthalene[p8[ . 
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